Generic Hubbard model description of semiconductor quantum dot spin qubits 
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We introduce a Hubbard model as the simple quantum generalization of the classical capacitance 
circuit model to study semiconductor quantum-dot spin qubits. We prove theoretically that our 
model is equivalent to the usual capacitance circuit model in the absence of quantum fluctuations. 
However, our model naturally includes quantum effects such as hopping and spin exchange. The 
parameters of the generalized Hubbard model can either be directly read off from the experimental 
plot of the stability diagram or be calculated from the microscopic theory, establishing a quantitative 
connection between the two. We show that, while the main topology of the charge stability diagram 
is determined by the ratio between inter-site and on-site Coulomb repulsion, fine details of the 
stability diagram reveal information about quantum effects. Extracting quantum information from 
experiments using our Hubbard model approach is simple, but would require the measurement 
resolution to increase by an order of magnitude. 

PACS numbers: 73.21. La, 03.67.Lx, 71.10.-w, 73.23.Hk 
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Introduction- Electron spin qubits in coupled quantum- 
dot systems are among the most promising candidates 
for quantum information processingi"~ because of their 
long coherence times^ and easy scalabilityi^ Two physical 
effects dominate the electronic properties of quantum- 
dot systems: the Coulomb repulsion and the quantum 
fluctuations associated with electron hopping and spin 
exchange ji^—iS 

Due to Coulomb repulsion, the number of electrons on 
each dot can be controlled by gate voltagesj^ This is 
mapped out by the charge stability diagram, which il- 
lustrates the equilibrium charge configurations as a func- 
tion of the two gate voltages on the dots. A classical 
capacitance circuit model neglecting all quantum effects, 
but including the Coulomb energy accurately, is almost 
universally used to model the charge stability diagram in 
quantum-dot systems. 

Quantum effects are essential for encoding and ma- 
nipulating quantum information in quantum-dot sys- 
tems. In the simplest scheme, which makes use of single- 
electron spin-up (-down) states, the exchange interac- 
tion constitutes two-qubit operations.^ In the double- 
dot system, the singlet (triplet) states form the qubit, 
the level splitting of which is controlled by the exchange 
interactiouf^"^ 

Since the charge stability diagram invariably serves as 
the experimental starting point for subsequent qubit ma- 
nipulations, understanding it in depth is of great impor- 
tance. Although the classical Coulomb effect is crucial 
in controlling the charge stability diagram, the quantum 
fluctuation effect should also manifest itself. A generic 
theory must be able to reconcile both. Previously, the 
charge stability diagram has been studied using the clas- 
sical capacitance model,— where the Coulomb inter- 
action is parametrized by the electrostatic energy of ef- 
fective capacitors. Although this method indeed pro- 
duces the main experimental features, ^'^ the quantum ef- 
fects are completely neglected in the capacitance model, 
which severely limits the applicability of the method in 



understanding quantum aspects of the charge stability 
diagram. 

In this Rapid Communication, we introduce a Hubbard 
model, as the simplest quantum generalization of the 
classical capacitance model, including all possible quan- 
tum effects allowed by symmetry, and develop a quantum 
theory for the charge stability diagram. The Hubbard 
model has earlier been discussed in the literatur o^^i^^ 
in the context of quantum-dot physics, but a complete 
quantitative description of the charge stability diagram 
connecting to the capacitance circuit model and the mi- 
croscopic confinement potential, as we provide here, has 
not been done before. We show that, in the absence 
of quantum effects, the generalized Hubbard model is 
equivalent to the capacitance model. However, while the 
quantum effects are intrinsically not allowed in the ca- 
pacitance model, they can be straightforwardly accom- 
modated in the generalized Hubbard model, making it 
suitable as a starting point of a more general theory in- 
cluding quantum fluctuations. We also note that some 
of the quantum effects were previously considered using 
a quantum mechanical two-level model)^ which can be 
naturally derived from our more general Hamiltonian by 
projecting onto the single-electron subspace. 

Microscopically, the quantum effects in the generalized 
Hubbard model originate from the overlap of the electron 
wave functions in a given confinement potential, which is 
treated using the configuration interaction method-^"^ 
For the bi-quadratic confinement potential)^ we iden- 
tify an important dimensionless parameter r] as the ra- 
tio between the height of the potential barrier and half 
the harmonic-oscillator energy-level spacing, which com- 
pletely determines the parameters of the Hubbard model 
as well as the shape of the charge stability diagram. Since 
these parameters can be extracted from a comparison 
with the experimental results, the generalized Hubbard 
model provides a quantitative bridge between the exper- 
iment and the microscopic theory. 

Mapping of capacitance model to Hubbard model- We 
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consider a double quantum-dot system with each dot (la- 
beled by z = 1,2) capable of holding iV^ = 0, 1 and 2 
electrons. Dots 1 and 2 are connected to the left and 
right gates with voltages Vl and Vr, respectively. 

We start with the classical capacitance model/ ^ in 
which the two dots are connected with a capacitance Cm 
and each dot is connected to the nearest gate via ca- 
pacitance Ci^/fj. The capacitance model assumes exact 
integer number of electrons residing on each dot, which 
we denote by {Ni,N2). The important quantity is the 
electrostatic energy as a function of Ni and N^^: 



E(Ni,N2) 



- EviNi ~Ev2N2 + Eq. 



Cm /C; 

{Cr^l + Cm) /Cl- Moreover, Eyia 
■2 r<2 



Here, the charging energies are Ec 

|2 , ^ ^ 



and Eci 2 

\e\{CLCRVL^R + CLCmVL+CnCmVR)/Cl, and C| = 
ChCm + CRCm + ClCr. Note that a precise value of the 
overall constant energy shift Eq is unimportant for the 
study of charge stability diagram. 

We map this model to the extended Hubbard model^- 
(without hopping), defined as 



H 



E 

1=1,2 



{-iliiii + UiUi^nii) + Ui2nin2, (2) 



where 



c\^Cia is the number operator of electron 



with spin a on site i, Ui — Uif + riii, iii and Ui denote 
the chemical potential and Coulomb interaction on site i, 
respectively. U12 denotes inter-site Coulomb interaction 
between dots 1 and 2. Since ni and n2 are good quantum 
numbers, the eigenenergy can be expressed as 



E' (ni,n2) 



1=1,2 



-fJ'int + —ni{ni 



1) 



+ Ui2nin2, (3) 



where the on-site Coulomb repulsion is present only 
when the dot is doubly occupied. {pl,^J-2) are linear 
combinations of {Vl,Vr) plus a constant energy shift^"^: 
Ail = |e| {aiVL + I3iVr) + 71, fi2 = \e\ {P2VL + a2^i?) + 
72. With this correspondence, comparing Eqs. ([T]) and 
Eq. ([3]), a mapping is immediately found: Ui = Eci, 

U12 = Ecm, ai,2 = {ClCr + CLM^m) /Cl^, "fi = 

—Ui/2. In particular, when Ui ^ U2 = U, we have 



1 - /3i,2 - ai,2 = U/{U + U12) = a. 



(4) 



Calculation of the charge stability diagram from 
Eq. ([2]) is straightforward. For given /ii and /i2, one 
identifies the ground state [labeled by electron occupancy 
(ni,n2)] by finding the lowest eigenenergy. The geome- 
try of the stability diagram is completely determined by 
U and U12, therefore one can read them off directly from 
the experimental plot. First, the length of the phase 
boundary between (1,1) and (0,2) is ■\/2C/i2, which de- 
fines Ui2- Second, the slope of the phase boundary be- 
tween (1,1) and (1,2) is a/[a — 1) and then U is ob- 
tained from Eq. (jj]). Figure [Ua) illustrates this idea. 
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FIG. 1: (Color online) (a) The charge stability diagram of 
the double dot system calculated using the extended Hub- 
bard model [Eq. 0], using parameters extracted from the 
experimental plot* which has been overlaid in the lower right 
part of the figure. U = 6.1 meV and U12 = 2.5 meV. (b) 
Charge stability diagram calculated with the same U and U12 
but with the hopping t varied, plotted in the /ii-/i2 plane. 
t = 0: black sohd fines, t = 0.6 meV: red dotted lines. Note 
that the main efi'ect of the quantum fluctuation is a rounding 
of the sharp triple point structure of the stability diagram. 



showing the charge stability diagram calculated from the 
extended Hubbard model, using parameters directly ex- 
tracted from the experiment^ which has been overlaid 
on the figure. The honeycomb shape is a direct conse- 
quence of a finite U12, and therefore the extended Hub- 
bard model [Eq. ([2])] is the minimal model that fits the 
experiment. 

Generalized Hubbard model- Quantum fluctuations ne- 
cessitate going beyond the extended Hubbard model. A 
general Hamiltonian for the double-dot system including 
up to two-body interactions can be written as 



H 



Fkiclci + Ghjkicic\ckCi (5) 

hjkl 



where the subscripts include site and spin indices. In this 
paper, we focus on the case without magnetic field. All 
possible terms that conserve the total particle number 
N — ni -\- n2 and the total spin Sz can be classified into 
the following categories: 

H = H^ + Ht+Hu + Hj (6) 

where the chemical potential part is H^^ = Y^^^i—^J■inia), 
and the hopping terms are Ht = J^ai^'^^ia^^a + H.c). 
The Coulomb repulsion terms are Hu = UiUi^nii + 
U2n2tn2i + C/i2("-it"24, + "-u"-2t) + {U12 - Je)ini-fn2^ + 
nnn2i), and Hj = - JecJ^4t'^2;Cii- - Jpcl^cl^cifCii - 

J2icr Jtniac\-^C2a + ii-C- , including the spin-exchange ( Je), 
pair- hopping (Jp), and occupation- modulated hopping 
terms (Jt)^^- In the following, we assume Ui = U2 = U . 
Equation ([6]) defines our proposed generalized Hubbard 
model. 

When Ht and Hj are present, ni and 712 are no 
longer good quantum numbers; therefore, the capaci- 
tance model ceases to be valid. The Hamiltonian is, 
in general, a 16 x 16 matrix for the double-dot system 
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but can be block diagonalized according to {N, Sz}- For 
states with different N, the phase boundaries are defined 
as before. For states within the same N block, the ground 
state is a superposition of (ni, 712) states and we label a 
phase according to the one that dominates the ground 
state. 

Figure mb) shows the calculated charge stability dia- 
gram with two different values of t (with = Jp = Jt = 
0), plotted on the ^1-/^2 plane. The t = data are the 
same as that shown in Fig.[T](a). Finite t changes the cur- 
vature of various phase boundaries and thus smoothens 
the sharp corner of the tripl-point structure of the stabil- 
ity diagram. For instance, the phase separator between 
the (0,0) and (1,0)/(0,1) complex has been changed to 
hyperbola /ii/i2 = t^- Because of the smallness of t and 
the insufficient resolution of the experiment in Fig. [IJa) , 
the value of t is not readily available from experiments 
at this stage, but t can, in principle, be read off from the 
details of the Vl-Vr diagrams. 

The generalized Hubbard model has several tunable 
parameters. However, we have a caveat: the parameters 
can not be arbitrarily chosen, otherwise unphysical sit- 
uations may occur. For example, if C/12 > U, the (1,1) 
state is completely ruled out from the stability diagram. 
Therefore, the model must be backed up by a microscopic 
theory leading to a reasonable choice of the parameter set 
defining the Hubbard model. 

Microscopic theory- In the microscopic theory^^ii^ one 
typically solves the many-electron Schrodinger equation 
with the confinement potential V (r). The potential wells 
near the centers of the quantum dots are approximated 
as parabolic, so the low-lying excited eigenstates of an 
isolated single dot are of harmonic-oscillator type, i.e., 
the Fock-Darwin states. The Fock-Darwin state on dot i 
can be written as \ipi^ni,mi {fi))^ with the principal quan- 
tum number rii = 0, 1, • • • and the azimuthal quantum 
number ruj = —Uj, —Uj + 2, ■ ■ ■ ^rij — 2^nj, the meaning 
of which is analogous to that of atomic orbitals. 

There are an infinite number of orbitals, but, for 
simplicity, we keep the lowest one denoted as S or- 
bital. The situation with multiple orbitals will be dis- 
cussed elsewhere. The Fock-Darwin states for differ- 
ent dots are, in general, not orthogonal; therefore, we 
build a new set of orthogonal basis by making the 
transformation {{'^i) {^2))^ = O^^^^ {\ipi) \(p2))^ , where 
O is the overlap matrix (Oki — {^k\ '-Pi)) generated 
by the Fock-Darwin states in the single particle sub- 
space. We recognize that the new basis \^i) actu- 
ally corresponds to c\ |0). The coupling parameters in 
Eq. ® are Fu = / dr-^l{r)h{r)-^ i{r), and Gum = 
/ dri J dr2'i/lir i)^* {r2)Hc'i'k{ri)'i'i{r2) where h{r) = 

[p — eA (r)]^ / {2m*) + V (r) is the single-particle Hamil- 
tonian, and He = /{Aneeolri — r2\) is the Coulomb 
interaction, with the effective mass of electrons m* = 
0.067me and relative dielectric constant e = 13.1 in GaAs 
quantum dot systems. 

In this calculation, we use a bi-quadratic poten- 
tial shown in Fig. V{x,y) = Min[l^(-o,0) - 
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FIG. 2: (Color online) (a) Schematic illustration for the bi- 
quadratic potential profile. 2a denotes the inter-dot distance, 
while vt represents the height of the central potential barrier 
when Hi = fJ.2 = /i. (b)-(f) Parameters in the generalized 
Hubbard model as functions of 77 = mujoa^/h, calculated at 
a — 15nm (black solid lines), 25nm (red dashed lines) and 
35nm (blue dotted lines). In panels (d) and (e), the calculated 
parameters do not change with a so only the black solid lines 
are shown. 



fii,W{a,0) — /i2,0], where the parabolic potential cen- 
tered at (xo,yo) is W{xo,yo) = muj^[{x - xq)'^ + (y - 
?/o)^]/2. The important parameters are the inter-dot dis- 
tance 2a and a dimensionless parameter r] = mujQa'^/h, 
denoting the ratio of the height of the central barrier 
(wf, = muj^a^ 12) and half the harmonic oscillator energy 
level spacing {huJo/2) when /ii = ^2- 

Figure [5] shows the microscopically calculated param- 
eters of the generalized Hubbard model as functions of rj 
for three different values of a. We note that the 77 < 1 case 
should be treated with caution: The central barrier is too 
low relative to the zero-point energy, effectively merging 
the two dots into one with very strong quantum fluctua- 
tions. Hereafter, we focus on the 77 > 1 regime. The on- 
site Coulomb interaction U is relevant for the extent of 
the localization of electron wave functions, controlled by 
ujQ, while the inter-site Coulomb interaction U12 depends 
on the inter-dot distance 2a. Figure HKb) and[5Jc) show 
these dependences: For a given a, increasing r] means 
that the electron wave functions become more localized, 
thus, U increases correspondingly, while U12 changes very 
slightly. If one fixes 77 but increases a (meaning ujq must 
be decreased), both U and U12 decrease, as expected. 

Unlike U and U12, which are constants over the entire 
stability diagram, t slightly changes as the potential is 
varied. Thus we take an average value t over the range 
that we consider: Vl,Vr £ [0,{U + 2Ui2 + 2fiwo)/|e|] 
(see Fig. 131). t should decrease either as the dots become 
further apart from each other, or as the electron wave 
function in each dot becomes more localized with in- 
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FIG. 3: (Color online) (a)-(c) Charge stability diagram cal- 
culated at three different values of r\ with a fixed as 33.8 nm. 
(d) The probability of electron configurations (f , f ) and (0, 2) 
along the blue dashed lines on panels (a)-(c) as a function 
of the detuning energy e = Vr — Vl- Black solid lines: the 
probability of (1,1) state; red dotted lines; the probability of 
(0,2) state. 

creasing. Figure [2{f) shows the ratio of t and /7 as a 
function of rj, and one can see that for rj > l,t decreases 
as either r/ or a is increased. Figures [H^d) andl^Je) show 
Ui2, Je{= Jp), and Jt divided by U, respectively. The 
calculations have been carried out for all three values 
of a but, remarkably, the results do not depend on any 
particular choice of a, leaving only a dependence on the 
dimensionless microscopic parameter rj. 

As mentioned above, U12/U plays a key role in deter- 
mining the shape of the stability diagram, with Je and 
Jt being an order of magnitude smaller than Ui2. This 
can be seen in Fig. [3l which plots the stability diagrams 
for three different values of 77. If the central potential 
barrier is low [Fig. ^a)] the two dots effectively merge 
into one, and the curvature of the phase boundaries in 
the neighborhood of the triple point is very small. For an 
intermediate value of ij [Fig.[31[b)], the stability diagram 



has the honeycomb pattern. As 1] further increases, the 
central barrier is high such that the overlap between the 
two dots is suppressed and the system behaves as two 
independent dots, leading to a tilted checker board pat- 
tern [Fig. [3l^c)]. Therefore, our approach gives a quanti- 
tative description of the evolution of the charge stability 
diagram as the inter-dot quantum fluctuation is varied, 
which can be compared to experiments^^ upon appropri- 
ately choosing parameters As noted above, the bound- 
aries between different configurations with the same N 
are determined by the probability of configurations. As 
an example. Fig. ^d) shows the probabilities of (1,1) 
and (0,2) states along the lines indicated as dashed lines 
in Figures |31[a)-[3Uc), as functions of the detuning energy 
defined a.s e = Vr — Vl- The crossing point defines the 
phase separator, and the width of the crossover encodes 
the quantum fluctuations. When r/ is small, the inter-dot 
coupling is large: the states are substantially mixed and 
the crossover is rather smooth. If 77 is large, the inter- 
dot coupling is small: the crossover is much sharper, as 
expected. Therefore, an alternative strategy of extract- 
ing quantum effects from experiments is to measure the 
width of the probability crossover near the phase separa- 
tor rather than from the curvature of the phase boundary. 
Our technique can be applied quantitatively to experi- 
mental results with sufficient resolution. 

Conclusion- We present a generalized Hubbard model 
approach to the quantum-dot system to study the charge 
stability diagram including quantum fluctuations. We 
identify the capacitance model as the limiting case of 
the generalized Hubbard model with vanishing quantum 
fluctuations. We show that quantum effects have observ- 
able consequences on the charge stability diagram, which 
can be extracted using the Hubbard model. We estab- 
lish a direct correspondence between the charge stability 
diagram and the microscopic double-dot confinement po- 
tential through the generalized Hubbard model, and show 
that the finer details of the stability diagram provide in- 
formation about the underlying quantum effects. 
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